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■ Based on the previously suggested model of nanoscale dislocations induced Josephson junctions 
and their arrays, we study the magnetic field induced electric polarization effects in intrinsically 
granular superconductors. In addition to a new phenomenon of chemomagnetoelectricity, the model 
predicts also a few other interesting effects, including charge analogues of Meissner paramagnetism 

' (at low fields) and "fishtail" anomaly (at high fields). The conditions under which these effects can 

be experimentally measured in non-stoichiometric high- Tc superconductors are discussed. 
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. Both granular superconductors and artificially prepared arrays of Josephson junctions (JJAs) proved useful in 
"j^ I studying the numerous quantum (charging) effects in these interesting systems, including blockade of Cooper pair 
tunneling [1], Bloch oscillations [2], propagation of quantum ballistic vortices [3], spin-tunneling related effects using 
specially designed SFS*- type junctions [4,5], novel Coulomb effects in SINIS-type nanoscale junctions [6], and recently 
. observed geometric quantization phenomena [7] (see, e.g., Ref. [8] for the recent review on charge and spin effects in 
' mesoscopic 2D Josephson junctions). 

More recently, it was realized that JJAs can be also used as quantum channels to transfer quantum information 
between distant sites [9-11] through the implementation of the so-called superconducting qubits which take advantage 
of both charge and phase degrees of freedom (see, e.g., Ref. [12] for a review on quantum-state engineering with 
^ Josephson-j unction devices). 

At the same time, imaging of the granular structure in underdoped Bi2Sr2CaCu20s+s crystals [13], revealed an 
\^ apparent charge segregation of its electronic structure into superconducting domains (of the order of a few nanometers) 
^-H , located in an electronically distinct background. In particular, it was found that at low levels of hole doping {5 < 0.2), 
the holes become concentrated at certain hole-rich domains. Tunneling between such domains leads to intrinsic 
. granular superconductivity (GS) in high-Tc superconductors (HTS). As was shown earlier [14], GS based phenomena 
' can shed some light on the origin and evolution of the so-called paramagnetic Meissner effect (PME) which manifests 
"p^ . itself both in high-Tc and conventional superconductors [15,16]. 

■ In this paper, within a previously suggested [14] model of JJAs created by a regular 2D network of twin-boundary 
^ , (TB) dislocations with strain fields acting as an insulating barrier between hole-rich domains in underdoped crystals, 

I , ' we address yet another class of interesting phenomena which are actually dual to the chemomagnetic effects described in 
^ , Ref. [14]. Specifically, we discuss a possible existence of a non-zero electric polarization P{B, S) (chemomagnetoelectic 
Q ■ effect) and the related change of the charge balance in intrinsically granular non-stoichiometric material under the 
Q influence of an applied magnetic field. In particular, we predict an anomalous low-field magnetic behavior of the 
effective junction charge Q{B,6) and concomitant magnetocapacitance C{B,S) in paramagnetic Meissner phase and 
a charge analog of "fishtail" -like anomaly at high magnetic fields. 



II. THE MODEL 



Recall that the observed [13,17-20] in HTS single crystals regular 2D dislocation networks of oxygen depleted 
regions with the size do of a few Burgers vectors can provide quite a realistic possibility for existence of 2D Josephson 
network within CuO plane [21,22]. In this regard, it is also important to mention the pioneer works by Khaikin and 
Khlyustikov [23-25] on twinning-induced superconductivity in dislocated crystals. 

At the same time, in underdoped crystals there is a real possibility to facilitate oxygen transport via the so-called 
osmotic mechanism [14,19,20,26] which relates a local value of the chemical potential /i(x) = pt(0) + V/x • x with a 
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local concentration of point defects as follows c(x) = e^A'(x)/fei3T^ allows us to explicitly incorporate the oxygen 
deficiency parameter 6 into our model by relating it to the excess oxygen concentration of vacancies c„ = c(0) as follows 
S = 1 — c„. Assuming the following connection between the variation of mechanical and chemical properties of planar 
defects, namely /i(x) = KClQe{x) (where e(x screened strain field produced by tetragonal regions in d- 

wave orthorhombic YBCO crystal, Q.q is an effective atomic volume of the vacancy, and K is the bulk elastic modulus), 
we can study the properties of TB induced JJs under intrinsic chemical pressure V/U (created by the variation of the 
oxygen doping parameter 5). More specifically, a single SIS type junction (comprising a Josephson network) is formed 
around TB due to a local depression of the superconducting order parameter A(x) oc e(x) over distance do producing 
thus a weak link with Josephson coupling J{5) = e(x)Jo = Jo{6)e~^^^/'^° where Jo{5) = cqJo = {Hv/ KQ,q)Jo (here 
Jo oc Aq/ Rn with Rn being a resistance of the junction). Notice that, in accordance with the observations, for 
stoichiometric situation (when 5 ~ 0), the Josephson coupling J (6) ~ and the system loses its explicitly granular 
signature. 

To describe the influence of chemomagnetic effects on charge balance of an intrinsically granular superconductor, 
we employ a model of 2D overdamped Josephson junction array which is based on the well known Hamiltonian 

n^f2Jij{i-cosct>ij) + f2^ (1) 

ij ij 

and introduces a short-range (nearest-neighbor) interaction between N junctions (which are formed around oxygen- 
rich superconducting areas with phases 0,), arranged in a two-dimensional (2D) lattice with coordinates Xi = (xj, t/j). 
The areas are separated by oxygen-poor insulating boundaries (created by TB strain fields e(xy )) producing a short- 
range Josephson coupling = Jo(<5)e^l''" 1/*^. Thus, typically for granular superconductors, the Josephson energy of 
the array varies exponentially with the distance Xy = Xj — Xj between neighboring junctions (with d being an average 
junction size). As usual, the second term in the rhs of Eq.(l) accounts for Coulomb effects where qi = —2eni is the 
junction charge with rii being the pair number operator. Naturally, the same strain fields e(xjj) will be responsible 
for dielectric properties of oxygen-depleted regions as well via the ^-dependent capacitance tensor Cij{S) = C[e(xjj)]. 

If, in addition to the chemical pressure Vy^(x) — iirr2oVe(x), the network of superconducting grains is under the 
influence of an applied frustrating magnetic field B, the total phase difference through the contact reads 

<^ij = 4>% + ^ (xii A Hi, ) • B + ^ ^ , (2) 

where ^9- is the initial phase difference (see below), rij^ = Xy/j | with Xjj = (xj -|- Xj)/2, and w = 2Xl{T) + 1 
with Xl being the London penetration depth of superconducting area and I an insulator thickness (which, within the 
discussed here scenario, is simply equal to the TB thickness [26]). 

As usual, to safely neglect the influence of the self-field effects in a real material, the corresponding Josephson 
penetration length Xj = i/$o/27r/Ioi^ must be larger than the junction size d. Here jc is the critical current density 
of superconducting (hole-rich) area. As we shall see below, this condition is rather well satisfied for HTS single 
crystals. 



III. CHEMOMAGNETOELECTRICITY 



In what follows, we shall be interested in the behavior of magnetic field induced electric polarization (chemomagne- 
toelectricity) in chemically induced GS described by a 2D JJA. Recall that a conventional (zero-field) pair polarization 
operator within the model under discussion reads [27,28] 

N 
i=l 

In view of Eqs.(l)-(3), and taking into account a usual "phase- number" commutation relation, [^i,nj] = i6ij, it can 
be shown that the evolution of the pair polarization operator is determined via the equation of motion 

dp 1 2e ^ 

= ^ [P> ^] = ^ Jij sin <P^j {t)'^ij (4) 

ij 
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Resolving the above equation, we arrive at the following net value of the magnetic-field induced longitudinal (along 
a;-axis) electric polarization P(5, B) =< Px{t) > and the corresponding effective junction charge 



T Z 

Q(5,B) = J dt J dt' J ^sm^{^,t')xe-\^\/'', 





(5) 



where S = 27rd^ is properly defined normalization area, r is a characteristic time (see Discussion), and we made a 
usual substitution Ylij ^ij(t) ^ I c'^2;^(x, t) valid in the long- wavelength approximation [28]. 

To capture the very essence of the superconducting analog of the chemomagnetoelectric effect, in what follows we 
assume for simplicity that a stoichiometric sample (with 5 ~ 0) does not possess any spontaneous polarization at 
zero magnetic field, that is P(0,0) = 0. According to Eq.(5), this condition implies = 27rm for the initial phase 
difference with m = 0, ±1, ±2, ... 

Taking the applied magnetic field along the c-axis (and normal to the CuO plane), that is B = (0, 0, B), we obtain 
finally 

g(j,5) = go(<5)^^±^M^ (6) 

^ ^ ^ ^(l + 62)(l + 62)2 ^ ' 

for the magnetic field behavior of the effective junction charge in chemically induced granular superconductors. 

Here Qo{S) = eTJo{5)/h with Jq{5) defined earlier, h = B/Bq, b = b — b^, and 5^ = B^/Bq ~ {ksTT/TijS where 
Bfj_{S) = {fiyT/fi)Bo is the chemically-induced contribution (which disappears in optimally doped systems with ^ ~ 0), 
and Bo = ^o/wd is a characteristic Josephson field. 

Fig.l shows changes of the initial (stoichiometric) effective junction charge Q (solid line) with oxygen deficiency 
5. Notice a sign change of Q (dotted and dashed lines) driven by non-zero values of 6 at low magnetic fields (a 
charge analog of chemically induced PME). According to Eq.(6), the effective charge changes its sign as soon as the 
chemomagnetic contribution B^((5) exceeds an applied magnetic field B (see Discussion). 

At the same time. Fig. 2 presents a true chemoelectric effect with concentration (deficiency) induced effective junction 
charge Q{S, 0) in zero magnetic field. Notice that Q{6, 0) exhibits a maximum around 5c — 0.2 (in agreement with the 
classical percolative behavior observed in non-stoichiometric YBa2Cu30r-s samples [17]). 

It is of interest also to consider the magnetic field behavior of the concomitant effective flux capacitance C = 
TdQ{d,B)/d^ which in view of Eq.(6) reads 

where $ = SB, and Ca{d) = tQo('5)/$o- 

Fig. 3 depicts the behavior of the effective flux capacitance C{S,B) in applied magnetic field for different values of 
oxygen deflciency parameter: 5 ~ (solid line), 6 — 0.1 (dashed line), and S = 0.2 (dotted line). Notice a decrease of 
magnetocapacitance amplitude and its peak shifting with increase of 6 and sign change at low magnetic fields which 
is another manifestation of the charge analog of chemically induced PME (Cf. Fig.l). 



IV. CHARGE ANALOG OF "FISHTAIL" ANOMALY 



So far, we neglected a possible field dependence of the chemical potential of oxygen vacancies. Recall, however, 
that in high enough applied magnetic fields B, the field-induced change of the chemical potential A^v{B) = ^v{B) — 
liv{0) becomes tangible and should be taken into account [14,29,30]. As a result, we end up with a superconducting 
analog of the so-called magnetoconcentration effect [14] with field induced creation of oxygen vacancies Cv{B) = 
Cy{Q)exp{—/S.iJLy[B)/kBT) which in turn brings about a "fishtail" -like behavior of the high-field chemomagnetization 
(see Ref. [14] for more details). 

Fig. 4 shows the field behavior of the effective junction charge in the presence of the above-mentioned magnetocon- 
centration effect. As it is clearly seen, Q{6{B), B) exhibits a "fishtail" -like anomaly typical for previously discussed [14] 
chemomagnetization in underdoped crystals with intragrain granularity (for symmetry and better visual effect we also 
plotted —Q{5{B),B) in the same figure). This more complex structure of the effective charge appears when the 
applied magnetic field B matches an intrinsic chemomagnetic field B^{S{B)) (which now also depends on B via the 
magnetoconcentration effect). Notice that a "fishtail" structure of Q{S{B),B) manifests itself even at zero values of 
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field-free deficiency parameter (5(0) (soUd line in Fig. 4) thus confirming a field-induced nature of intrinsic granular- 
ity [13,17-20]. Likewise, Fig. 5 depicts the evolution of the eflfective flux capacitance C{6{B),B) in applied magnetic 
field B/Bq in the presence of magnetoconcentration effect (Cf. Fig.3). 



V. DISCUSSION 



Thus, the present model predicts appearance of two interrelated phenomena (dual to the previously discussed 
behavior of chemomagnetizm [14]), namely a charge analog of Meissner paramagnetism at low fields and a charge 
analog of "fishtail" anomaly at high fields. To see whether these effects can be actually observed in a real material, 
let us estimate an order of magnitude of the main model parameters. 

Using typical [17,19] for HTS single crystals values of — 150nm, d ~ lOnm, and jc — 10^'^^/m^, we arrive 

at the following estimates of the characteristic Bq — 0.5T and chemomagnetic B^{S) ~ 0.5Bq fields, respectively. So, 
the predicted charge analog of PME should be observable for applied magnetic fields B < 0.25T. Notice that, for the 
above set of parameters, the Josephson length is of the order of Aj ~ Ifim, which means that the assumed in this 
paper small-junction approximation is valid and the "self-field" effects can be safely neglected. 

Furthermore, the characteristic frequencies w ~ r"^ needed to probe the suggested here effects are related to the 
processes governed by tunneling relaxation times r ~ %/Jq{5). Since for oxygen deficiency parameter 6 = 0.1 the 
chemically-induced zero-temperature Josephson energy in non-stoichiometric YBCO single crystals is of the order 
of Jo ((5) — ksTcS ~ ImeV, we arrive at the required frequencies of u; ^ 10^^ Hz and at the following estimates 
of the effective junction charge Qo — e = 1.6 x 10~^^C and flux capacitance Cq ~ 10~^*F. Notice that the above 
estimates fall into the range of parameters used in typical experiments for studying the single-electron tunneling effects 
both in JJs and JJAs [1,2,12,31] suggesting thus quite an optimistic possibility to observe the above-predicted field 
induced effects experimentally in non-stoichiometric superconductors with pronounced networks of planar defects or 
in artificially prepared JJAs. (It is worth mentioning that a somewhat similar behavior of the magnetic field induced 
charge and related flux capacitance has been observed in 2D electron systems [32].) 

And finally, it can be easily verified that, in view of Eqs.(l)-(5), the field-induced Coulomb energy of the oxygen- 
depleted region within our model is given by 



with Q{6,B) and C{6,B) defined by Eqs. (6) and (7), respectively. 

A thorough analysis of the above expression reveals that in the PME state (when B <C B^^) the chemically- 
induced granular superconductor is always in the so-called Coulomb blockade regime (with Ec > Jo), while in 
the "fishtail" state (for B > B^) the energy balance tips in favor of tunneling (with Eq < Jo)- In particular, 
Ec{5^ B = Q.IB^) = -I Jo ((5) and Ec{5, B = B^) = ^Jo{5). It woifld be also interesting to check this phenomenon of 
field-induced weakening of the Coulomb blockade experimentally. 



In conclusion, within a realistic model of 2D Josephson junction arrays created by 2D network of twin boundary 
dislocations (with strain fields acting as an insulating barrier between hole-rich domains in underdoped crystals), 
a few novel electric polarization related effects expected to occur in intrinsically granular material under applied 
magnetic fields were predicted, including a phenomenon of chemomagnetoclcctricity, an anomalous low-field magnetic 
behavior of the effective junction charge (and flux capacitance) in paramagnetic Meissner phase and a charge analog 
of "fishtail" -like anomaly at high magnetic fields as well as field-dependent weakening of the chemically-induced 
Coulomb blockade. The experimental conditions needed to observe the predicted here effects in non-stoichiometric 
high-Tc superconductors were discussed. 
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FIG. 1. The effective junction charge Q{5, B) /Q{5c,0) (chemomagnetoelectric effect) as a function of apphed magnetic field 
B/Bo, according to Eq.(6), for different values of oxygen deficiency parameter: (5 ~ (solid line), S = 0.1 (dashed line), and 
S = 0.2 (dotted line). 
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FIG. 2. Chemically induced effective junction charge Q(S,0)/Q{Sc,0) in a zero applied magnetic field (true chemoelectric 
effect). 
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FIG. 4. A "fishtail" -like behavior of effective charge Q{S{B), B)/Q{5c,0) in apphed magnetic field B/Bq in the presence 
of magnetoconccntration effect (with field-induced oxygen vacancies S(B)) for three values of field-free deficiency parameter: 
5(0) ~ (solid line), 5(0) = 0.1 (dashed line), and 5(0) = 0.2 (dotted line). 
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FIG. 5. The behavior of the effective flux capacitance C{S{B) , B) /C{Sc, 0) in apphed magnetic field B/Bo in the presence of 
magnetoconcentration effect for three values of field- free deficiency parameter: 5{0) ~ (solid line), 5{0) = 0.1 (dashed line), 
and 5(0) = 0.2 (dotted line). 
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